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Abstract 

The Greer, Pujo-Menjouet and Webb model [Greer et at, J. Thcorct. Biol., 242 (2006), 598-606] 
for prion dynamics was found to be in good agreement with experimental observations under 
no-flow conditions. The objective of this work is to generalize the problem to the framework of 
general polymerization-fragmentation under flow motion, motivated by the fact that laboratory 
work often involves prion dynamics under flow conditions in order to observe faster processes. 
Moreover, understanding and modelling the microstructurc influence of macroscopically moni- 
tored non-Newtonian behaviour is crucial for sensor design, with the goal to provide practical 
information about ongoing molecular evolution. This paper's results can then be considered 
as one step in the mathematical understanding of such models, namely the proof of positivity 
and existence of solutions in suitable functional spaces. To that purpose, we introduce a new 
model based on the rigid-rod polymer theory to account for the polymer dynamics under flow 
conditions. As expected, when applied to the prion problem, in the absence of motion it reduces 
to that in Greer et al. (2006). At the heart of any polymer kinetical theory there is a configura- 
tional probability diffusion partial differential equation (PDE) of Fokker-Planck-Smoluchowski 
type. The main mathematical result of this paper is the proof of existence of positive solutions 
to the aforementioned PDE for a class of flows of practical interest, taking into account the flow 
induced splitting/lengthening of polymers in general, and prions in particular. 
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1 Introduction 



1.1 Taking space into account for our problem: what is new in biology, 
what is new in mathematics? 

In 1999, Mascl et al. [12] introduced a new model of polymerization in order to quantify some 
kinetic parameters of prion replication. This work was based on a deterministic discrete model 
developed into an infinite system of ordinary differential equations, one for each possible fibril 
length. In 2006, Greer et al. in [6] modified this model to create a continuum of possible 
fibril lengths described by a partial differential equation coupled with an ordinary differential 
equation. This approach appeared to be "conceptually more accessible and mathematically 
more tractable with only six parameters, each of which having a biological interpretation" [6]. 
However, based on discussions with biologists, it appeared that these models were not well 
adapted for in vitro experiments. In these experiments, proteins are put in tubes and shaken 
permanently throughout the experiment to induce an artificial splitting in order to accelerate the 
polymerization-fragmentation mechanism. To the best of our knowledge, dependence of polymer 
and monomer interaction on the shaking orientation and strength, space competition and fluid 
viscosity had never been taken into account until now. Thus, it seemed natural to propose a 
model generalizing the Greer model and adapt it to the specific expectations of the biologists. 

We therefore introduce a new model of polymer and monomer interacting in a fluid, with the 
whole system subjected to motion. A large range of in vitro experiments involving this protein 
refers to this protocol in order to accelerate the polymerization-fragmentation process. Moreover, 
even as our model could be well adapted to other polymer-monomer interaction studies, we give 
here a specific application to prion dynamics to make an interesting link with the previous Mascl 
et al. [12] and Greer et al. [6] models. On the other hand, due to the complexity of the model, any 
mathematical analysis becomes a challenge. We adapt here a technique of semi-discretization in 
time for proving the main result of existence of positive solutions, we also provide the basis for the 
numerical approximation of the problem. The mathematical novelty of this paper resides in the 
choice of the ad hoc function spaces and the appropiate modification of the existing techniques 
to this new type of problem. Also this work presents an alternative way for proving the existence 
of positive solutions as compared to the one given by Engler et al. in [5], Laurengot and Walker 
in [11] and Simonett and Walker in [17]. It is then useful to those who consider which techniques 
to use when proving the existence of positive solutions of this class of equations. 

The objective of this paper is twofold: not only to make a step forward in mathematical 
modelling of a class of polymer-monomer interaction models, but also to propose, within a new 
framework, how to adapt an existing mathematical technique that will prove the existence of 
positive solutions to the problem. The biological implications {e.g. quantitative and qualitative 
comparison with experimental data) of this paper model will be addressed in a subsequent work. 
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1.2 The polymer-monomer interaction model: an application to prion 
dynamics 

Prion proliferation is challenging at both the biological and mathematical levels. Prions are re- 
sponsible for several diseases such as bovine spongiform encephalopathy, Creutzf eld- Jacob disease. 
Kuru and it is now commonly accepted that prions are proteins [14]. 

For the sake of clarity, we present several fundamental morphological features of prions with 
relevance to the mathematical modelling of this paper {i.e. molecular dynamics of a low enough 
concentration prion solution) . 

There are two types of prions: the Prion Protein Cellular also called PrP c and Prion Protein 
Scrapie denoted by PrP Sc . It has been proven that PrP c proteins are naturally synthesized by 
mammalian cells and consist only of monomers. On the other hand, the infectious PrP Sc proteins 
are present only in pathologically altered cells and exist only in "polymer" -shape. The conversion 
process of a non-pathological into a pathologically modified one consists in attaching the former to 
an already existing polymer (for details see e.g. [10]). As a consequence, the polymers lengthen. 
However the sized-up new polymers are fragile, and shorten down their size by splitting whenever 
the polymer solution is subjected to some flow conditions. The size lengthening/shortening 
process takes place continuously, its kinetics being dependent on monomer concentration, flow 
intensity, polymer size, etc. 

Polymers may be seen as string-like molecules [16]. When polymer proliferation occurs, they 
do interact to form fibrils; these latter exhibit a (physically speaking) more stable structure and 
appear as rod-like molecules (see figure 1). In this paper we deal with idealized rod- like PrP Sc , a 
realistic choice taking into account the flow-related experiments we investigate. We consider the 
presence of a finite amount of PrP c (free monomers) and PrP Sc proteins, as well as of "seeding" 
rod- like PrP at initial condition, and fibril lengthening/splitting (i.e. fragmentation). It is 
also important to note that our model is related to in vitro experiments: neither source terms of 
monomers and polymers nor degradation rates are taken into account. 

We propose a comprehensive molecular model that accounts for the flow behavior as observed 
in in vitro experiments, focusing on the dynamics of monomers and fibrils. A good deal of 
experimental laboratory work involves complex flows (e.g. diffusion, mixing, etc.). Raw data 
are provided by sensors designed to acquire macroscopically observable properties like stresses, 
flow rates, etc. The latter can strongly be influenced by the microscopic interactions. Our 
model docs provide an understanding of how various polymers-monomer and polymer-solvent 
relationship result in a configurational probability diffusion equation, with the help of which one 
can investigate the stress tensor and related quantities. Therefore, it is of use for flow pattern 
monitoring sensors. 

The current approach is at an early stage of development. The scission (breakage) process - 
the most important mechanism in the in vitro development/proliferation of infectious proteins 
- is taken three-dimensionally. While prior models such as those of [6, 12] (for mathematically 
in nature aspects related to, see [2, 5, 7, 11, 15, 16, 17, 18]) neglect the flow influence on prion 
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Figure 1: View of prion fibrils, Transmission Electron Microscopy image (Courtesy of Prof. J. -P. 
Liautard, inserm Universite Montpcllicr 2, France). 

dynamics, the one in [6] was rather succesful in predicting prion molecular dynamics in the in 
vivo rest state, and our model is a generalization of [6]. 

The prion fiber is modelled as a rigid rod polymer molecule the length of which is time 
dependent; sec figure 2. 

The dynamics of rigid rod molecular fluids has been initiated by Kirkwood [9] and signifi- 
cantly enriched and brought to fruition by Bird and his school [1] (see also [8] for a more succinct 
presentation) . As in any kinetical theory, the cornerstone is the probability of the configurational 
diffusion equation, which is of a Fokker-Planck-Smoluchowski type. The latter is the key ingre- 
dient for calculating (the macroscopic) stress tensor and related quantities. In the following we 
shall derive a suitable generalization of equations 14.2-8 in [1] that account for prion dynamics 
as observed in experiments [3, 14, 19]. 

This paper begins by first presenting the constitutive assumptions which later lead to the 
probability configurational diffusion equation in its general form. We give a mathematical con- 
ceptual framework and a presentation of the main result: the existence of global weak non- 
negative solution. To achieve this, we obtain a variational formulation of the corresponding 
boundary value problem, and the proof is based on a semi-discretization in time technique. The 
uniqueness of the solution will be proved in a subsequent paper. 
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Monomers 

Figure 2: Prion fibril modelized as rigid rod polymer under flow. 

1.3 The general model 
1.3.1 Polymers 

Let a fiber be modeled as a rod- like molecule here represented by a vector in R 3 . For convenience, 
we use separate symbols for the length r <E R + = (0, +oo), and for the angle- vector r\ £ 8 2 , with 
§ 2 being the unit sphere of R 3 . Contrary to the assumption made in [6] and for simplicity, 
we assume here that polymers could be arbitrary small, that is no critical (lower) length is 
considered (this assumption is explained in [4]). For technical reasons and without any loss of 
realistic assumptions, we suppose that fibers are contained in a bounded, smooth open set fl 
in R 3 , and the position of each fiber center of mass is denoted by the vector y. We assume a 
velocity vector field u : 57 x R + — > R 3 such that 



V y • u = in f2, and u • n = on <90. (1) 

with n the outward normal. The polymer configurational probability distribution function 
i/>(r, 77, y, t), at any time t > 0, solves the following equation 

d d 

— i> + u • V y + ^ (r(0, u, r, r))iP) = Btp + F^. (2) 

with (r, 77, y) £ R + x § 2 x fL Fibers are transported by the velocity vector field u and lengthening 
occurs at a rate r > that depends on the free monomers density, <f>. In dilute regime, the 
microscopic hydrodynamics is accounted for by the term B as in [13] and defined by 

[B$\(r, r/, y, t) = A(r) V„ • [DiV„V - (V y u r?) V] , (3) 
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where V r; and (V,,-) denote the gradient and divergence on § 2 . A > is a weight function 
that accounts for the influence of the length increase upon the motion and D\ > the diffusion 
coefficient on the sphere. Moreover, the transport on the sphere due to the velocity field is given 
by P ?) -l (V y u rf), with P ?) -lz = z — (z • T])i], for all z £ R 3 , denoting the projection of the vector 
z on the tangent space at rj. 

The fragmentation (scission) process takes place at rate /3(V y u, u, r, 77) > and is described by 
J- following [6] and given by 



00 



[Ftl>](r,r),y,t) = -pi) + 2 /?(V y u, u, r', r,) K (r, r')^(r' , r,, y, t) dr' . (4) 

J r 

The size redistribution kernel k accounts for the fact that a polymer breaks into smaller 
fibers. It is symmetric, since a polymer of size r' breaks with equal probability into a fiber of 
size r' — r and r; moreover, the fragmentation/recombination is mass preserving process. We 
assume here that upon splitting, given the pecularity of the motion process, and its impact on 
the scission, the resulting clusters of fibrils have the same center of mass as the initial polymer. 
It seems reasonable to assume that the orientation remains unchanged right after the scission. 
Therefore: n(r, r') > 0, /t(r, r') = if r > r' , n(r' — r, r') = «(r, r') and 

K(r,r') dr = 1. (5) 

The probability configurational function ifi must be a non-negative solution, satisfying the non- 
zero size boundary condition 

tp(o,v,y,t) = 0, (6) 

and the initial condition 

ip(r,i],y,0) = tp°{r,r],y), (7) 
with ip° a known non-negative initial probability. 

1.3.2 Monomers 

The concentration of free monomers, given by the distribution <p(y, t) at time t > at any y S Q, 
solves 



d 



u V y 0-L> 2 A0 = - / T{<f>,u,r,ri)ip(r,ri,y,t) dr dr], (8) 



with Z?2 > the diffusion coefficient. The integral term is due to polymerization of monomers, 
being transconformed (misfoldcd), into fibers. Moreover, monomer concentration <f> must be a 
non-negative solution satisfying the (no transport across) boundary condition 

V y 4> ■ ft = on dCl, (9) 



G 



with n the outward normal vector on the boundary dtt, as well as the initial condition 



0(y,O) = o (y), (10) 

with (fp an initially non-negative given concentration. We adjoin to these equations the balance 
equation for the total number of monomers contained in the domain O: 



<t>(y, t) + I r ip(r, t], y, t) drj dr 

,+ xS 2 



dy = p, for all t > 0, (11) 



where p is (experimentally) known from the outset. The above balance equation is formally 
satisfied, as a consequence of equations (2)-(8) using also (1). 

1.3.3 Velocity vector field and momentum balance equations 

As an aside, notice the velocity vector field, u(t, y) S R 3 , for all t > and y G il, satisfies the 
Navier-Stokes equations (for incompressible fluids) 

d 

— u + (u • V) u = - Vp + fAu - V • S, 

V-u = 0, (12) 
u • ft = 0. 

p is the pressure , v the viscosity of the Newtonian solvent within which the prions (i.e. rigid- rod 
molecules) are dissolved, and S is the non-Newtonian extra stress tensor contribution (to the 
total stress) due to the presence of rigid rods. The latter is given by [1] as 

S(y,i) = r 2 rj®rj Tp drjdr. (13) 
Jr + Js 2 

In this paper, we suppose that u is given and the unknown functions are only ip and <f>. The 
existence and uniqueness of the solutions to the full system with the Navier-Stokes equations 
introduced above (that is u, ip and <f) will be the topic of a subsequent paper. 

1.4 Constitutive assumptions 

Assume the velocity vector field satisfies the regularity 

ueC 1 ([0,oo),F^ 1 ' oo (r!)) (14) 

such that 

V y • u = and u • n = on dfl. (15) 
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Next, we adhere to the view on prion proliferation expressed in [6, 7, 12, 15]. The splitting 
(scission) rate of fibers, given by /3, is assumed to be linear in r. Therefore let g : M3(R) xl 3 x 
§ 2 — > R + be continuous with respect to the first and second variable, such that (3(a, v, r, rf) = 
g(a,v,ri) r, for all a £ M^(M.), v £ R 3 , r > and rj £ § 2 . Moreover, we assume that for all 
bounded subsets B C R 3 and O C Ma(R) there exist positive constants ~g B > g g such that 

g B O < g((T, v, rj) < g B Ol for every (a, v, 77) £ O x B x § 2 . (16) 

Let T > be fixed. Then, due to the smoothness of u, there exists ~g > g > such that 

for every (t,y, V ) £ [0,T] xSlxS 2 , g < g(V y u, u, rj) < g. (17) 

We consider the polymerization rate r linear in (the free monomers density) <p, i.e. there exists 
To > such that 

T(<j),v,r,r)) =T <j). (18) 

This assumption had been already evoked by Greer et al. [6] and corresponds to a mass action 
binding. The splitting kernel k accounts for the probability of a polymer with initial length r, 
to split into a polymer with a shorter length r' as described in [6], and is given by 

fl/r' ifO<r<r', 
K {r,r') = { ' (19) 
[0 else. 

This expression is compatible with (5) (and the conservation law (11)). Then the length weight 
function A > is supposed to be in L°°(R + ) and there exists Ca > such that 

HIU-(O) = C A < 00 (20) 

We remark that, by virtue of u being sufficiently smooth and for fixed T > 0, there exists Cp > 
such that 

l|FV ( v y u V) l|L°°([o,T]xnxS 2 ) = C P < 00, (21) 
Using the result stated in the Appendix, there exists Cd > such that 

||V„ • P^-l (V y u ?y) || L «([o,T]xnxS2) = C D < 00. (22) 
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Thanks to the assumptions given in this section, the problem can be re-written as: 
d d 

— V> + u • V y ^ + To^—ip - A(r) V n • [£>iV^ - P v ± (V y u r?) V] 

/•DO 

= -.9(V y u, u, 17)7^ + 2. 9 (V y u, u, 77) / ^(r' ; 77, y, i) dr', 

St 



(23a) 



u • V y (f> - D 2 A<t> = -T (f> f ip(r,r),y,t) dr dr), (23b) 

JS 2 xR + 

ip(r = 0,r],y,t) = 0, (23c) 

V.y0 • rl = 0, on dil (23d) 

ip(t = 0) = V>° and 0(t = 0) = 0°, (23e) 



1.5 Particular case: zero velocity field, as in the Greer's model 

Consider u = 0, and assume that g is such that g(0,r)) = go, a constant, for any rj. In fact, 
even in the absence of flow the prion-fibrils can undergo scission and re-combination. Suppose 
that 4> is independent of y, then let f(t,r) = —— Jq xS 2 ^(^i^jY)*) drjdy be the average of ip. 
Integrating equations (23) leads to 



g- f f + rom^f + 9orf = 2gj /(/, t) dr' over (t, r) e 

±-6{t) = -T*4>{t) I f(r,t)dr, (24) 
at J M+ 

l/(0,*) = 0. 

Note that the above system of equations is the one proposed in [6] where it was produced under 
the assumption of prion conservation mass (no protein synthesis, no metabolic degradation). 

2 Variational formulation and main result 

First we present the functional framework one of the main mathematical novelty of this paper, 
next the definition of weak solutions to the system (23) , and eventually the proof of the existence 
of a weak solution of this system. 



9 



2.1 Functional framework 



Let a : M + ->■ M + be defined by a(r) = c al for a a > 0. Denote Q = § 2 x : 
Let the following Hilbcrt spaces be defined as 



Then, 



and 



L 2 a = UeL 1 loc (nxQ), f ^dqdyKoo) 

I JQxQ ) 



and dq = a(r)drdrj. 



(25) 



V = \i> e iL (0 x Q) , / (A(r)|V^f + (1 + r)^ 2 ) dqdy < oo 

JaxQ 



V x = UeL] 0C {£lxQ), 



d \ 

—ip + A{r)\V v ^\ 2 + (1 + r)i/> 2 dqdy < oo 

Recall the Sobolev space ff 1 (J7) endowed with the norm 



(26) 



(27) 



HWm = \\4>\\z»(a) + \\Vv4>\\v>w- (28) 

We also use the canonical embedding 

V 1 C V C L 2 a = (L 2 J cV'c {Vx)' . (29) 

For any 6 <E K, let = G L i 1 oc (!) X Q) , /q x q M r e drdrjdy < ooj. Then we have the 
canonical embedding 

L 2 C Lj, for any a > and > 0, (30) 

which makes sense in regard to the mass conservation and the total quantity of polymers when 
= or = 1. 



2.2 Variational formulation 

To begin with, we introduce test function spaces. Let T > 0. First, for the polymer -0-equation, 
let X\ be the completion of C^°((— T, T) x Q x S 2 x [0, +oo)) with respect to the norm || • \\x 1 
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a 





L 2 



Ml 



dt 



(31) 



In particular, this implies that, if tp 6 X\, then ip(t = T) = 0. Second, the test functions for the 
0-equation are elements of X2, the latter space being the completion of C^°((— T, T) x f2) with 
respect to the norm 7J 1 ((0, T) x f2). In particular this implies that if <fi € X2, then <j>(t = T) = 0. 
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In order to obtain a variational formulation of (23) we first assume that we have a strong solution 
which is smooth enough. Then we multiply (23a) by ip(r, f], y, t)a(r), with ip € X\, and integrate 
over (0,T) x Q, x Q, next we multiply (23b) by <\> € X2 and integrate over (0,T) x 0. We note 



To</'^'0 i> a(r)dr = - j r # ^ (^ a M) 



dr. 



( J-V 1 + cr0 



+ 



(32) 



since ip E Xi. One also has: 

V„ • (Z>iV„V) $dr] = - J DiV„V • (v„$ - 2#) (33) 



s 2 



and 



/ V„ • (P^ (V y u rj) i>) i) dn = - f P n ± (V y u 77) $ ■ h v i> - 277-0) drj, 

JS? 7§2 \ / 

= - / *V (V y u 77) ■(/> • d?7, 
Js 2 



(34) 



since P„-l (V y u 77) ■ r/ = (see for instance Appendix II in [13] for calculation details on the 
sphere). Moreover, by assumption (15) on u, 



(u-V y ^dy = - / V(u-V y V) dy, (35) 

J n 

and 

f (u-V y( />)0dy = - / 0(u-V y 0) dy. (36) 



Then a variational formulation of (23a) is 

- / 1>° $(t = 0) dqdy - [ f r/> + u • V y /) dqdy dt 

Jqxq Jo Jnxo / 



JflxQ 
T 



A(r) (DiWr,^ (Vriip - 2r/rpj - P v ± (V y u 77) i< ■ V v rpj dqdy dt 
I I i'(^(^y u > u 'V)r$ -T (j>(JI-^+a^^J dqdy dt 

= 2 f [ ,g(V y u, u, r?) ( [ V dA V~ dqdy dt, 

JO JQxQ \Jr / 



for any ip £ X\, 
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and for (23b), 



0° 4>(t = 0)dy- 
n Jo Jn 

T 

D 2 V y • V y + tqc. 



<p + u - V y </> I dy dt 



o Jn 



ip drdr) 



S 2 xl 



dy dt = 0, 



(38) 



for any <p E X 2 . 

2.3 Main result: existence of non-negative solutions of the problem 

At this point we are prepared to introduce our main result. It gives the existence of non- negative 
weak solution to our problem under the general assumptions of section 1.4. 

Theorem 2.1 (Main result). Let 0° <E be non-negative and ip° £ L 2 a non-negative such 

that there exists a constant Co > with 

i>° < C Q e~ ar . 



Then, for any T > 0, there exists at least one solution {ip, (f>) to the weak formulation (37)-(38) of 
the problem (23), with ip and <fi non-negative. Moreover we have ip € L°°(0, T; L^) n L 2 (0, T; V) 
and <P g L°°(0,T; L 2 (n)) H L 2 (0,T; H 1 ^)). 

Remark 1. Proving the uniqueness of the solution is a rather lengthy undertaking and will be 
done in a follow up paper. 

Remark 2. ; Weak solutions to the above variational formulation with stronger regularity than 
the one implied by the theorem above satisfy the problem (23) in a strong sense. Moreover, 
this variational formulation complies weakly with the mass conservation principle. Therefore, let 
if e F^O.T), with(p(t = T) = 0, and take ip{r,n,y,t) = re- ar ip(t) E X x and(f(t,y) = ip(t) £ X 2 
in the variational formulations. Using the fact that, for any real value function f 



[ V ■ V„/ dq = 0. 

J§ 2 



(39) 



we obtain 



Mt = o) 



+ / r iP° dr, dr 

JR + xS 2 
T j 



T Mt) 
o ai Jn 



r ip dr\ dr 



_xS 2 



(40) 



dy dt = 0. 



If the solution is smooth enough we have then the mass conservation result 



dt 



+ I r ip dr] dr 

l + xS 2 



dy = Q. 



(41) 
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3 Proof of the main result 



The proof consists of three main steps. First (subsection 3.1), a semi-discretization in time of 
the problem to obtain an approximation of the solution. Second, we get appropriate estimates 
(subsection 3.2), and third we obtain a solution by passing to the limit (subsection 3.3). 

3.1 Semi-discretization in time 

Let N > and {t n }^ =0 a subdivision of [0, T] such that to — 0, tpj = T and t n — t n -\ = At > 0. 
We denote by ip n and 4> n the approximations of ip and <j) at t n . Denote u"(y) = u(t n ,y). First, 
for any s E [0, T], consider the following problem on [0,T]: 

4-X™(*) = u n (x n (t)), 

dt A w yA y ,h (42) 

x(s) = y- 

We recall that the regularity of u is C^O, T; W 1 ^), therefore u" 6 W l '°°(Q.) so that there exists a 
unique solution x™ which will be denoted in the following by s , y)- The map y — > x"(i; s, y) 
is a homeomorphism from Q onto J7, and since u is divergence-free, we have 

dct V yX n (t; s, •) = 1, a.c. in f2 x [0, T]. (43) 

Define the function 

x n : n -> 0, by x„(y) = X n (i n ;*n-l,y)- (44) 
This map i„ is invertible. Let us denote z„ as its inverse. We remark that 

z«(y) = x n (i„_i;f„,y). (45) 

Assume now that ^i™ -1 6 1/ and £ -ff 1 arc known. We consider two problems: 
find ip n E y such that 

j> n (r,V,y) ~ ^ n ~ 1 (r,y,z n (y)) ~ 



(46) 



-if 9 



2 / 5 (V y u", u", r?) f /°° ^ dr') $ a(r)drd V dy, 

JflxQ \Jr / 
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for any ip G V\, and find <j) n £ H 1 such that 
> "(y)-0"- 1 (y) 

D 2 V y 0" • V y + T O( /> 



for any 4> € H 1 . Problem (46) is re- written as 

a"(V",T/0 = W), for any ^ 6 Vi 



with 



a" = a 1 " + a 2n 



where a , a arc defined on V x Vi by 



a ln (^i,<<92) = / 4(r)DiV,^i ■ iy v <P2 - ^Wz) dqdy 
JnxQ 

- / A(r)P n ± (V y u n 7 ? ) v?i • V,,^2 rfqrfy 



TQ 



( J^¥>2 + OL ip 2 



+ / .g(V y u n ,u",?7)r^i^2 dqdy 



and 



At 



(fHP2 dqdy, 



nxQ 



respectively, and is defined on l? a by 



2 f 5 (V y u", u", V ) ( ^ dA if dqdy 

JClxQ \Jr J 



— I ip n 1 o z n ip dqdy. 



At 



nxQ 



The problem (47) is re-written as 



for any 0eH\ 
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with b n defined on H 1 x H 1 such that 

b n (<p U (p 2 ) 



J n ( Al^ 1 ^ + ' VyV?1 ) V2 + ° 2 Vy(pi ' y y^ 2 ) dy 



V> n_1 drdri I dy, 



(54) 



and I? defined on I? by 



(55) 



Lemma 3.1. Let N £ N*, </>° £ </>° > and ip° £ L 2 a such that 

< V° < C e- ar a.ein Q 

with Cq > a constant. 

Then there exist two sequences {i[> n }n=i C and {0"}^Li C iJ 1 (Q) satisfying (48) and (53). 
Moreover, for At small enough, we have that: 



and 



and 



< V" < Cooe" 
O<0" < ll0 O ||r 



/or every n £ {0, 1, • • • iV}, 
for every n £ {0, 1, • • ■ iV}, 



max 

n=0,—N 



[ \r\ 2 dqdy + ^AiV / A(r)|V^"| 2 dqdy 
+2gAi V / r\^ n \ 2 dqdy + V / |V> n - ^Z-"" 1 ° z n \ 2 dqdy 

n=1 JflxQ n=1 JQxQ 

<ie k ^U°f Li 



max 

n=0,---N 



' N N 

/ |0"| 2 dy + V / \ ( j> n -<f> n - 1 \ 2 dy + 2D 2 Aty2 / |V y 0"| 2 dy 
Jo „ =1 Jn n=1 Jn 



< 



2H \\h m , 



(56a) 
(56b) 



(57) 



(58) 
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where in the above we denoted 



ki = — , 
a 



k 2 = ar \\<j) Q +C d C a , 
Coo = 2C e^ +k ^ T , 



and 



(Recall Cd cind Ca are given by equations (20) and (22)). 



Proof of Lemma 3.1 

Let us consider the sequence of numbers {C n }^ =0 defined by induction as 

C' n = \ + r^ Cn-ii for every n = 1, ■ ■ ■ N. (59) 
with Co as in the hypothesis of the Lemma. 

We proceed by induction. Suppose that and <f> n ~ 1 are defined as elements of V and L°°(D,), 

respectively. Suppose also that 

< V™ -1 < C„_ie~ Qr , (60a) 
< < (60b) 

We shall prove the existence of ip n £ V and </>" 6 iJ 1 (Sl) solutions of (48) and (53), respectively. 
We also prove that they satisfy 

< V™ < C n e~ ar , (61a) 
<<j> n < (61b) 

The above inequalities give (56a) and (56b) since we have 

for At small enough. 

Step 1. Regularization and existence. 

We introduce a regularization of a", denoted a" defined on V\ x Vi, 

f d d 

a"(</5i,^ 2 ) = e / ^-</?i^-</?2 rfqrfy + a n (^i,^ 2 )- (63) 
JnxQ or or 
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We shall first prove the existence of a sequence (ip™) E in V\ solutions of 

a?(i%,$) = IZ$), forany^GFi (64) 

Clearly a™ is bilinear and continuous on V\ X V\. Next we prove the coercivity of a™. Indeed, let 
ip GVi and we remark that 



/ 2t] ■ V v p ip dii = i] ■ V^p 2 di] = 
Js 2 Js 2 

since ■ rj = 2 and r\ ■ r\ = 1. One has 



(65) 



f \A(r)P n ± (V y u t?) ip ■ V n <p\ drj < \ [ [B^lVr, 

Js, 2 z Js 2 \ 



v{ 2 + C^A^2 



Finally, 



+ 



to / </>"~ W</> a(r)dr < --ar / <j> n ~ l p 2 a(r)dr. (67) 



We remark that this inequality can be proved by using a regularized sequence {<Pm) m that 
converges to p in V\ and the fact that the remaining term in the right-hand side of (67) can 
be dropped according to its appropriate sign. Then, invoking (60b) and the above remarks, it 
follows that 

al n (p,p)>^- [ A(r)|V^| 2 dqdy + g [ rp 2 dqdy 
z JnxQ JnxQ 

- -L (aroD^WLoo + C 2 P C A ) [ p 2 dqdy, 
ZL> i JnxQ 

which in turn implies 



(68) 



c£(p>,p) > e 



Jnxi 



d 2 
to* 



D f 

dqdy+^- A{r)\V v p\ 2 dqdy 

1 JnxQ 



2Di V At 

The coercivity of a" follows for At small enough. 
Next, due to the inequality (60a), we have 



g rp 2 dqdy (69) 

JnxQ 

aro^ill^llioo - C 2 P C A / p 2 dqdy, 



inxc 



< (70) 



P CO 

/ ip^dr' 

■J r 

which implies that, for any cp € 

[ g(\7 y u n ,u n , V )( ^ dr'\ \p\ dqdy < [ \p\ drd V dy. (71) 
JnxQ \Jr J a JnxQ 
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One also obtains 

/ f^ozn M dqdy < C n _! f \ip\ drdrjdy. (72) 
JnxQ JnxQ 

We deduce that Z™ £ (La)' C (Vi)' by the continuous embedding of in L 1 . Applying the 
Lax-Milgram theorem, for all e > there exists a unique G V\ solution of (64). Next we will 
prove the existence of solutions to (53). First, b n is clearly a bilinear and continuous function on 
H 1 x H 1 . To prove its coercivity, let ip 6 H 1 . Since 

/ u" • S7 yV <p = I / U " • V yV ? 2 = (73) 



wc have 



>]- [ V 2 dy + D 2 f |V y ^| 2 dy, (74) 



using the positivity of tp n x , and thus b n is coercive. Moreover, l£ £ (-Z? 1 )' since <j> n 1 £ L°°. As 
a consequence of the Lax-Milgram theorem, there exists a unique 4> n G H 1 satifying (53). 

Step 2. L°° - Estimates 

To begin we first prove two estimates for ip™: for its y-norm and for its derivative with respect 
to r. It follows from (69) and the continuity of Z™ that there exists a constant C > 0, dependent 
of At, such that 

(A(r)|v,^"| 2 + (i + r)|vrf) < c, 

(75) 

dqdy < C. 



e 



9 it 



'QxQ 

Next we prove the non- negativity of -0™ and 4> n . Let us denote [•]+ and [•]_ respectively the 
positive and negative part, both positive valued. Then, <f> n = [<f> n ]+ — [4> n \- and these two parts 
belong to H 1 . We have 

iKm-) = b n (r, r]_) = -b n ([<t> n ]-, m-) m 

and invoking (55) and (60b), Z£([0 n ]_) > 0. Therefore 

b n (ir]-,W l ]-)<0, (77) 
hence <j? L > 0. Next, -0" = — ['0™]-' the positive and negative parts belong V\, and 

= <(c m\-) = -<mu m~), m 

Invoking (52) and (60a), Z"([V>™]-) > 0. Thus 

m~) < 0, (79) 
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hence ip™ > 0- Let us now obtain L°° estimates . We have, according to (60b) and using the 
above notation, that 

6 n ([^-||^ IU-]+ ) [^ B -ll^ lk-]+) 

= b n (r-n o \\L~d0 n -H o h~] + ) 

= b n (<t> n , W 1 - II0°IM+) - W l - U°\\l~] + ) (80 ) 

- W* - H0°lk~]+) - 6 n (ll^°IU», [0" - 

Then by (60b) 

^(r-ii^ o iu-] + ,[0 n -ii0 o iu-] + )<o, (si) 

hence 4> n < ||0 o ||i,°°- Let C„ as defined in (59); then 

aUm-C n e- ar } + , [r e - C n e- ar }+) 

= <(^ 1 - C n e- ar , - C n e- ar }+) 

= - C n e~ ar ] + ) - <(C„e- Qr , [V£ - C n e~ ar } + ) 

= ~ C n e~ ar } + ) ~ a?(C n e- ar , [V£ - C n e~ ar ] + ). 



Next, for any ip £ V\ positive, 

a"(C„e~ Qr , <p)=-e [ aC n -^-ip drdrjdy 
JQxQ ar 



(82) 



/ C n A(r)P v ± (V y u"r/) • V v <p drdrjdy 

JUxQ 

[ C n T 4> n ~ 1 ( -^-(f + a ip ) drdrjdy 
JnxQ \ dr J 

/ C n g(V y u n ,u n ,r])ripdrdT]dy + C n — I <p drdrjdy. 

JQxQ At JflxQ 



(83) 



We remark that 



; / aC n -^-ip drd-qdy = -e f aC n tp(r = 0,r?,y) drjdy < 0, 

JQxQ ar JQxS 2 

[ CnTor' 1 ^ drdridy = - f C^^cpir = 0, V ,y) d V dy < 0. 
JnxQ clr JnxS 2 



(84) 



19 



Then, by (20), (22), (60b) and the positivity of ip, 

a n e {C n e- a \ ip) > [ C n A(r)V v ■ (P^ (V y u" ?7 )) <p drdndy 



+ C n (J^ - ar ||</>°|| L ^ f drdndy 

^ C n \ T7 - fc 2 ) / f drdr]dy. 



(85) 



Moreover, by (52), (71) and (72) 

IM < C n -1 + V drdndy. (86) 

\a AtJ J nxQ 

Now, replacing ip by [ip™ — C„e~ Qr ] + and using (82) (85) and (86) one gets 



a n s ([^-C n e- ar } + , [C - C n e- ar }+) 



< 



C n _i ( fci + — | - C n f — - fc 2 



t/j drdrjdy. 



(87) 



At y "V At 

Using now the particular form of C n gives 

" C„e~ Q 1 + , [VC " C„e-H + ) < 0, 
hence 

W < C n e~ ar . (89) 

Step 3. Convergence and positivity 

The sequence (ipe) e obtained for all e > is uniformly bounded in V by (75), so it weakly 
converges to an element ip n G V up to a subsequence. Moreover, (e 1 ^ 2 ^:^?) s is bounded in i^, 
then for e — >• 0, i/i" solves (48). The positivity of ip™ yields the positivity of ijj n . Moreover, by 
virtue of (89), ip n for e — > 0, and inequalities (56a) are satisfied. 

Step 4. Additional estimates 

From (69), (52) and (56a) one gets 

= <(V £ n , VC) " We) >^ / A(r)\Vr,^\ 2 d q dy 

z JnxQ 

+ g r\^\ 2 dq.dy 

/ Qx r Q (90) 

™3 



/ ICI 2 dqdy 
JnxQ 
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Remarking that 2si(si — S2) = s\ + (si — s 2 ) 2 — s\ for any reals s±, S2, leads to 
Di / A(r)|V^ e "| 2 dqdy + 2ff / r|^"| 2 dqdy 



+ / Dd 2 + IC - V-™- 1 o z n \ 2 - ^ o z„| 2 ] dqdy (91) 
m -/nxQ 



< fc 3 / |^| 2 dqdy. 
JSlxQ 



Then, taking the liminf for e — > 0, multiplying by At and using the fact that 
/ |^- 1 oz„| 2 =/ IV.- 1 ! 2 , gives 



DiAt [ A(r)\\7 v i> n \ 2 dqdy + 2gAt[ r\ip n \ 2 dqdy 

JnxQ JnxQ 

+ (1 - k 3 At) f \i> n \ 2 dqdy + [ |V" - <F _1 o z n \ 2 dqdy (92) 
JnxQ JnxQ 

< f l^"" 1 ! 2 dqdy. 
JnxQ 

Multiply the last inequality by (1 — k^At) 11 " 1 and sum over n from n = 1 to n = N. Use the 
inequality 

(1 - k 3 At) n > (1 - fc 3 A<) w > X -e- k * T 
to get (57). Taking <\> — <\> n in (47) and using (56b) and (73) we obtain 

J- / (|^|2 + |^-0»-i|2_|^-i|2) dy + ^ 2 / |V y 0"| 2 dy<O (93) 
Summing over n from 1 to JV produces (58), which ends the proof. 

3.2 Construction of a solution 

We now define, for any N large enough, the following functions 

^N(;t) = t —^r(-) + t -^r~ 1 (-), *e [*»-!,*»] (94) 

and 

TP+(;t) = ^ n (-), ^(;t)=i> n -\-), *G (*„-!,*«] (95) 

for n = 1, • • -N. 

We shall use analogous notations for </>jv and ujy. Let tp € Xi, (j) 6 A" 2 , both be test functions 
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and set tp = f t " ^ tp dt and 4> — f t " 1 4>dt. It is clear that G V\ and <f> G iJ 1 (f2). Then 



" a n (i> n ,i>(;t))dt= %(i>(;t))dt, 

t n -l Jtn-1 

" b n (r,$(;t))dt= /*" lUi>{;t))dt. 



Adding these inequalities, we obtain, for any %j) G X\, 

i>N{r,r],y,t) - i/j^(r,r),z N (y,t),t) 



o Jnxo At 



i/j(r,r),y,t) dqdy 



r f t N <y,t) J~ N (y,t) ^ t) dydf + r r {u + ) dydf 

Jo Ja At Jo Jn 

+ [ [ D 2 V y 0+ • V y dydt + r f [ 4>%\ [ ip N drdiA 4> dydt 
Jo Jn Jo Jn Ws 2 xR + / 



= 0. 



(96) 



+ D 1 f f A(r)V n ilj+ • (v„$ - 2rp£) dqdy 
Ja JqxQ k ' 

- f f A{r)P v ± (V y u+ V ) yj+ ■ dqdy (97) 

Jo JilxQ 

+ J ^ff(V y u^,u^,77)r^-T </>^ (Jj^^+a^jj dqdy 

= 2 / ff(V y u+ , u+ , rj) ( / ipjj dr' ) $ dqdy, 

Jo JflxQ \Jr / 

where in the above, 

Xjv(y,t) = x n (y) and z N (y, t) = z n (y), for any t G (t n -i, t n ). (98) 
Proceeding likewise, for any <f> £ X 2 , 



(99) 



However, to evaluate the limit At — > 0, we need some additional convergence results about the 
approximations. First, let us define the maps, 



Ai#](y,t)=/ ip(r,r),y,t) drdr), 

(100) 

A 2 [iP\(r, V ,y,t) = / i)(r', V ,y,t) dr', for any -0 6 L 2 (0,T;L 2 a ). 



We have the following lemma: 

Lemma 3.2. Let G L°°(Q,), <j)° > and ip° G sucA i/iai 
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< V° < Coe~ ar a.ein Q 



with Co > a constant. For {iPn}n o-nd {V'ivljv' constructed by virtue of Lemma 3.1, there 
exists ip G L 2 (0,T;V) (~l L°°(0,T;L 2 ) ; positive, such that, for N — > +oo we have the following 
convergence, up to a subsequence of N: 



A^ N ] - Ai[V] 



.1/2 







* - ttfeaJWy in L°°(0, T; 

weakly in L 2 (0 , T '; L a ) , 
weakly in L 2 (0 , T '; Z, a ) , 
weafc/y in L 2 ((0,T) x 0), 
weakly in L 2 (0 , T '; L 2 a ) . 



(101) 
(102) 
(103) 
(104) 
(105) 



Proof It is clear from (57) that 



and 



We then deduce that 



From (57) one infers 



ip% is bounded in L 2 (0, T; V) 



is bounded in L°°(0,T; L 2 a ). 



ij) N o z N is bounded in L°°(0, T; L 2 ). 



ipff — tp N ° z n in the norm of L 2 (0, T; L 2 a ) 



(106) 

(107) 
(108) 

(109) 



Then there exists i>+ 6 L 2 (0, T; V) n L°°(0, T; L 2 a ) and V" G £°°(0, T; L 2 ) such that, up to a 
subsequence in N we have 



ijfjj — 1 V> + weakly in L 2 (0,T;V) 



^n^^ * -weakly in L°°{0,T; LI), 



and 



■0jv o zjv — * -weakly in L°°(0,T;L 2 ). 



(110) 
(111) 

(112) 
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On the other hand we have 



Xn(y) -y = x n (t n ;t n - 1 ,y) - x n (V-i;*n-i,y) 
= At u"( x n (€;*n-i,y)). 

This implies 

||xiv(y,*) -y|Loc(]o,r[xn) ^ At H u lli-(]o,T[xr2) 
Now, for any %p € C£°(Q x fix]0,T[), with the help of (114) and (107), we obtain 



J!ix 
T 



[^(n^y,*) - ^Ar(^??,zjv(y,i),0] ^{r,r],y,t) dqdydt 



lo Jnx 
< CAt\\u 



^ N ( r ^ 7 h y, t) ip{r, r), y, t) - ip(r, rj, xjv (y, t),t) 



dqdydt 



L°°([o,T]xn) llrllc 



(113) 



(114) 



(115) 



We deduce that ip N — ip N o — > in the sense of distributions T)'(Qx]0, T[). This leads to the 
conclusion that = ip~ , and we denote by ?p the common value or ^ — . Therefore (101), 
(102) and (103) are proved. Let now ip <E L 2 ((0,T) x 0) 



(AiVjv _ A iV^) <p(y,t) rf y rf i 

o Jn 

/ tj}J l ipe~ ar dqdydt / ■>pLpe~ ar dqdydt 

o JnxQ Jo JnxQ 

-)■ 0, as TV +oo 



(116) 



since <^e _ar € L 2 . Now, invoking (101), proves (104). Finally, let ip G L 2 (0,T;L 2 ) and with the 
help of (103) we get 



(A 2 ^+ - A 2 i/>) %p dqdydt 

o JnxQ 

T r „ i-T , 

rtp^tp dqdydt — j j rtjjip dqdydt 



JQxQ 

0, as iV -> +oo. 



(117) 



JQxQ 



Which proves (105). The positivity of ?/> follows from the positivity of ip n for any n. This ends 
the proof. 

We now focus on the convergence of the 0jy sequence. 
Lemma 3.3. Let cp" e 0° > and ip° G L 2 smc/i i/iai 
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< V° < Coe~ ar a.ein Q 



with Co > a constant. For {4>n}n and N ; constructed by virtue of Lemma 3.1, there 

exists <j> £ L 2 (0,T; H 1 ) n L°°(0,T;L 2 ) positive such that we have the following convergence, up 
to a subsequence of N : 



Vy^Ar V y <j 

4>%, <S>n ->■ 



weafcZy L 2 (0,T;L 2 ) 
strongly L 2 (0,T; L 2 (fl)) 



(118) 
(119) 



Proof From (58), we deduce that 

0+ is bounded in L 2 (0,T; H 1 ^)), 



and 



Since we have 



we deduce that 



and 



'N 



ON 



b% is bounded in L°°(0, T; L 2 ) 



is bounded in L 2 (^,T;ff 1 (n)) for any £e]0,T[. 
, t n — t t — t n -i _)_ 

0jv is bounded in L°° (0, T; L 2 ) 

is bounded in L 2 (5, T;H\n)) for any 5e]0,T[. 



(120) 
(121) 

(122) 



(123) 



(124) 



It follows there exists a <f> 6 L 2 (0, T; if 1 ) n£°°(0, T; L 2 ) such that (118) is satisfied. On the other 
hand, from the equality 



; >/v 



in J.n—1 



dt At 

and from (47) we deduce that for any <fi £ H 1 ^) we have 



on [i„_i, t n ] 



(125) 



dt 



dy = - / u+ ■ V y (t> N 4> dy - D 2 I V y <^ • V y dy 



./n 



7"0 / 



Q WS 2 x 



(126) 



VV drdrj </> dy 



Using (120) and (107), gives 



'■>N 



dt 



is bounded in L 2 {Q,T; {H x {££))' 



(127) 
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Then, up to a subsequence of N, we have 



<j) N ^<P strongly in L 1 (6, T;I/(0)), for any 8 <=]0,T[. (128) 
Let us now prove that 

(j) N -> 4> strongly in L 2 (0, T; L 2 {il)). (129) 
We fix £ > and we have for any 8 G ]0, T[: 

HN-(f>\\hfn)dt= \\<Pn - <P\\ 2 L 2 {n) dt.+ Un - (f>\\ 2 L 2 (n) dt 
o Jo Js 



(130) 



<2CS + I \\(t>N-<t>\\l* (Q) dt 



5 



where C is an upper bound for ||0Af||L=°(o,T;L 2 ) and ||^||l~(o,T;L 2 )- Now taking 8 = we obtain 
from (128) that for N large enough 



T 



£ 



\\<t>N-<f>\\h m dt<-, (131) 
which proves (129). From (58) one gets 



Using the fact that 



and 



leads to 



This ends the proof. 



b+-^^0 strongly in L 2 (0,T;L 2 (ft)). (132) 



^-^N = t -^t^N-^N) ( 133 ) 



~ 4>n = t -^ ± {4>% ~ (134) 



N -4>%-+0 strongly in L 2 (0, T; L 2 (ft)). (135) 



3.3 Final stage of the proof of the main result 



In the following we let N -> +oo in (97) and (99) with %p E C~((-T,T) x Cl x § 2 x [0, +oo)) 
and <j> G C£°((— T, T) x fi x § 2 x [0, +oo)), respectively. We now prove that ijj and </> given by 
Lemmas 3.2 and 3.3 satisfy the variational equalities (37) and (38), respectively. Since At is 



2G 



small enough, we have 



— T7 ^(r,V,y,t) dqdydt 

f I-, ,x 17. xjy(y, t), t) - %j){r, n,y,t- At) , . 

/ y> N {r,V,y,t) t- dqdydt (136) 

JfixQ ^ 



1 

At 



At 

/ ip°(r, V, y)ip{r, V,y,t- At) dqdydt. 

JiixQ 



Smoothness of ip entails 

rAt 



1 

At./ 



V-°(r, r?, y)V(r, r?, y, t - At) dqdydt -> / ^<ty(* = 0) dqdy, (137) 
nxQ JnxQ 



and 



^0, y, *) ~ ^(r, V,y,t- At) 9~ 2 2 
At ► ^W^.y.*) strongly in L (0, T; £ Q ). (138) 

We also have 

^(r,7?,x n (y),t) -V>(r,77,y,i) ~, , , , . , , 
— = V v ip(r,r),y + 6>i(x„(y) -y),t) • £;v, (139) 

with 0i e]0, 1[ and 

f* = (140) 

Since x n (y) = x"(t„-i, t„, y) we have 

x"(t„-i,t n ,y) -x"(t„,t n ,y) 



At 

-^(t„_i + 2 At,t n ,y) = -u n (x n (t„_i + 2 At,t„,y)), 

at 



with 2 €]0,1[. Then 

^(r,r?,x n (y),t) - ^(r,n,y,t) 



On the other hand, for any s E [t n — i,t n ] 



x n (s;t n ,y)-y = x n (s;t n ,y)-x n (*« ; *„ . y ) 
<9% 



(141) 



At (142) 
= -V y $(r, i h y + 0i(x n (y) - y), t) • u™( X "(t„-i + 2 At, t„, y)). 



t (t n + 6 3 (s-t n ),t n ,y)(s-t n ), (143) 
u"(x"(t„ + 3 (s - t n ),t n ,y))(s - t„), 
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with 3 e]0, 1[, then 

\x n (s;t n ,y) -y| < |s-i„| l|u|U~(f2x]o,T[)- (144) 



Then one deduces from (142) and (144): 

i>(r, r}, x N (y, t),t)- V?(r, 77, y, t) 



->-u(t,y)-V y V(r,»?,y,t), (145) 



At 

strongly in L 2 (0,T;L 2 a ). Next, from (136), (137), (138) and (145) one gets 

\ — tt i>{r, T], y, t) dqdydt 

JilxQ n - t 

->■ - [ [ tp (^-4> + u ■ Vyip) dqdydt - I = 0) c 

Jo JnxQ \ot J JnxQ 

Now, from the strong convergences 



(146) 



V y u+ -> V y u, (147) 
.g(V y u+,u+,?7) ->#(V y u, 11,77), (148) 

and the fact that 



(149) 



one easily calculates the limit in (97) and gets (37). Moreover, 



fo,t) dydt 
iP°^(t = 0) dy - [ [ (f>^-4> dydt. 



(150) 



m Jo Jn dt 

Calculating the limit in (99) easily leads to (38). 

4 Conclusions 

Understanding polymer dynamics under different experimental conditions is of importance for 
the laboratory biologists. In this work we studied the influence of an external velocity field on the 
polymer-fibrils fragmentation (scission) and lengthening process. To the best of our knowledge 
this type of study has never been taken into account in the mathematical modelling of this 
problem. And even if our approach is at its early stage of development, we managed to obtain a 
rather good generalization of the existing models using more realistic assumptions when adapted 
to the prion study. 

In this work, we generalized the corresponding Fokkcr-Planck-Smoluchowski partial differen- 
tial equation for rigid rods in order to account for the fragmentation/lengthening process adapted 
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for prion proliferation. Moreover, we have introduced a set of two equations on monomers and 
polymers with a known flow. We prove existence and positivity of weak solutions to the system 
with assumptions on the rates and distribution kernel. The proof is based on variational formu- 
lation, a semi-discretization in time, and we obtain estimations which allow us to pass to the 
limit. To achieve this, we introduced a suitable functional framework (see section 2.1). 

The matter of existence of solutions to the full system (i.e. considering the time dependence 
of monomers together with the Navier-Stokes equations given in section 1.3) will be adressed in 
a future work. 
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Appendix 

Let M £ A^3(ffi), r\ G § 2 , we shall compute in spherical coordinates according to the base 

(eg,e<p, e r ) 

V v ■ P v ±Mrj = V v ■ Mr) - V n ■ (Mr) • 77)77. 
Note that in spherical coordinates, rj = e r and for F a vector value function, 

V„ • F = d e Fg + ^-F e + -^-:d v F v + 2F r , 
smO suit/ 

with Fk = F ■ efc, for k = 9,tp,r. According to the derivative of the vector of the base, see 
Appendix II [13] and the fact that 

dkMe r ■ ej = Mdy& r ■ e.j + Me r ■ d^e^ 

assumed that F = Me r , then 



V, ; • Me r = Me 6 



Next, take F = (Me r ■ e r )e r , it is clear that 

Fg = (Me r ■ e r )(e r ■ eg) = 0, and F v = (Me r ■ e r )(e r ■ e v ) = 0, 

thus 

• (Me r ■ e r )e r = 2Me r ■ e r . 
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Finally, 

V,, • P v ± M 77 = Meg ■ eg + M e v ■ e v - 2Me r ■ e r . 
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